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Abstract 

We present a Monte-Carlo generator for a variety of Glauber-like models (the 
wounded-nucleon model, binary collisions model, mixed model, model with hot 
spots). These models describe the early stages of relativistic heavy-ion collisions, 
in particular the spatial distribution of the transverse energy deposition which ul- 
timately leads to production of particles from the interaction region. The original 
geometric distribution of sources in the transverse plane can be superimposed with 
a statistical distribution simulating the dispersion in the generated transverse en- 
ergy in each individual collision. The program generates inter alia the fixed- axes 
(standard) and variable-axes (participant) two-dimensional profiles of the density 
of sources in the transverse plane and their azimuthal Fourier components. These 
profiles can be used in further analysis of physical phenomena, such as the jet 
quenching, event-by-event hydrodynamics, or analysis of the elliptic flow and its 
fluctuations. Characteristics of the event (multiplicities, eccentricities, Fourier co- 
efficients, etc.) are stored in a ROOT file and can be analyzed off-line. In particular, 
event-by-event studies can be carried out in a simple way. A number of ROOT scripts 
is provided for that purpose. Supplied variants of the code can also be used for the 
proton- nucleus and deuteron- nucleus collisions. 
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Program summary 

Title of the program: GLISSANDO October 2007, ver. 1.5 

Catalog identifier: 
Program summary URL: 

http:/ /www.pu. kielce.pl/homepages/mryb/GLISSANDO/index. html 
Program obtainable from: 

http:/ /www. pu.kielce.pl/homepages/mryb/GLISSANDO/index. html 
Licensing provisions: none 

Computer: any computer with a C++ compiler and the ROOT environment [1], 
tested with Intel Core 2 Duo T5200, 1.6 GHz, 1 GB RAM 
Operating system under which the program has been tested: Linux Ubuntu 7.04- 
8.04 (gcc 4.1.3-4.2.3), Scientific Linux Linux CERN 3.08 (gcc 3.2.3), Windows 
XP with Cygwin (gcc 3.4.4 cygwin special), MacBook Pro OSX 10.5.2 (gcc 
4.0.1), ROOT ver. 5.08-5.18 

Programming language used: C++ with the ROOT libraries 
Memory required to execute with typical data: below 120 MB 
No. of lines in distributed program, including test data: 2000 
No. of bytes in distributed program, including test data: 35 kB 
Distribution format: tar.gz 

Nature of physical problem: Glauber-like models of the initial state in rela- 
tivistic heavy-ion collisions 

Method of solution: Glauber Monte-Carlo simulation of collision events, an- 
alyzed with ROOT 

Restrictions concerning the complexity of the problem: none 
Typical running time: 40 s/10000 events for the wounded-nucleon model and 
60 s/10000 events for the mixed model with the T distribution, minimum-bias 
Au+Au collisions with the fixed-last algorithm for nuclear repulsion, disper- 
sion of sources, and hard-sphere wounding profile. A typical "physics" run with 
500000 events takes about 1 hour. The use of the Gaussian wounding profile in- 
creases the time about 6 times. (All times for Intel Core 2 Duo T5200 1.6 GHz) 



1 Introduction 



The Glauber Monte-Carlo simulations have become a basic tool in the analysis 
of relativistic heavy-ion collisions (for a review and the discussion of physics 
issues see the recent review [2] and references therein) . The approach provides 
the initial condition arising just after the collision of two relativistic nuclei. 
Within the semiclassical Glauber model, during the first stage of the collision, 
individual interactions between the nucleons deposit transverse energy. These 
elementary processes are classified as wounded nucleons or binary collisions. In 
this paper the individual deposition of the transverse energy at a definite po- 
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sition in the transverse plane is called a source and the distribution of all the 
sources created in an event, weighted by their individual deposited strength, is 
called the fireball. A weight, called relative deposited strength or RDS, 
is assigned to each source when the distribution of sources in the transverse 
plane is convoluted with a statistical distribution. This convolution simulates 
the dispersion in the generated transverse energy (particle production) from 
each elementary collision. The normalization of the deposited weight can be 
treated as a additional parameter. This parameter can be separated in the 
form of an overall factor multiplying the fireball density profile. In the follow- 
ing the deposited weight in each interaction is defined relatively to the one 
corresponding to an elementary nucleon-nucleon (NN) interaction. The RDS 
strength from an elementary interaction encompasses the contribution of two 
wounded nucleons and one collision and gives on average 1 (see section|3]). 

The Glauber Monte-Carlo methods are used to determine the centrality classes 
in the experiment, since the total nucleus-nucleus cross section is not measured 
at RHIC. Many analysis use the Glauber-model initial state, which can serve 
as input for cascade simulations or hydrodynamics. Moreover, the initial eccen- 
tricity of the shape of the fireball carries over to the elliptic flow coefficient v 2 
in the momenta of produced particles. The behavior of v 2 and its fluctuations 
is a sensitive probe of the dynamics. 

Although numerous codes exist for the Glauber Monte-Carlo calculations and 
the simplicity of the algorithm allows anybody to produce his own with no 
difficulty, we have decided to make ours publicly available for several rea- 
sons. Firstly, to our knowledge there is no published publicly available pack- 
age specifically dedicated to Glauber Monte-Carlo calculations. Most popular 
codes are much broader in their scope [3,4], having the initial-state generation 
only as one of the stages. This makes modifications, extensions, or compar- 
isons more difficult. Many details of the model and algorithm implementations 
are different in these codes, such as the choice of the nuclear density profiles, 
the variant of the Glauber model used, the implementation of nucleon-nucleon 
repulsion, etc. GLISSANDO offers flexibility here: 

• There is practically no constraint on the form of the nuclear density profile. 
In Sec. 12. II we discuss the subtleties related to the choice of the distribution 
of the centers of nucleons in the nucleus such that the profiles determined 
in the electron-nucleus scattering are reproduced. 

• We allow for superimposing a distribution of weights over the distribution of 
individual sources (wounded nucleons, binary collisions), since elementary 
collisions may result in the deposition of a varying amount of transverse 
energy. 

• Hard-sphere or Gaussian wounding profiles may be chosen (see Sect. 12.21) . 

• The position of the source may fluctuate relative to the center of the wounded 
nucleon or binary collision from which it originates, which accounts for the 
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finite size of the nucleons. 

• We offer two different algorithms to simulate the short-range nucleon-nucleon 
repulsion in the generation of the distribution of nuclei in the nucleon. 

In short, GLISSANDO tries to encompass all features and details of the used 
approaches. Moreover, 

• The variable- geometry analysis [5-8], accounting for the fluctuation of the 
center of mass and the direction of the axes of the moment of inertia of 
the fireball, is built-in in the code. This phenomenon is responsible for an 
increased eccentricity of the initial condition and thus is important in the 
studies of the elliptic flow. 

• We compute the two-dimensional density profiles for both the fixed-axes 
(standard) and variable-axes (participant) geometry, as well as the corre- 
sponding Fourier radial profiles. These are stored and can be used "off-line" 
in other physics analysis, such as the jet quenching, the initial condition for 
hydrodynamics, etc. 

• GLISSANDO can also be helpful for studies of event-by-event fluctuations of 
multiplicities, as measured for instance at the CERN SPS by the NA49 
collaboration [9]. 

• Basic characteristics of the events are also stored, allowing for event-by- 
event analysis of various quantities. 

• The code can also be used for proton-nucleus and deuteron-nucleus colli- 
sions. The Hulthen distribution is used to describe the NN separation in 
the deuteron. 

• The code is written in C++ and uses the ROOT libraries and data structures. 
That way it is tailored for the experimental community. We have taken 
an effort to document the program such that it can easily be modified or 
extended. 

• As an option, GLISSANDO generates the full event tree containing the po- 
sitions and RDS of the sources within fireballs. This information may be 
further processed "off-line" by other existing codes. 

We hope GLISSANDO will prove useful for heavy-ion physicists. 



2 Basic method 

The simulation of an event consists of three stages: 1) generation of the posi- 
tions of nucleons in the two nuclei according to the nuclear density distribution, 
2) generation the transverse positions of sources and their relative deposited 
strength (RDS) - see Sect. I3.1[ and 3) calculation of the event-by-event av- 
eraged physical quantities and fireball profiles and writing the results to the 
output file. Stage 3) is performed on-the-fly with stage 2). 
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2. 1 Generation of positions of nucleons 



The radial nuclear density distribution determined from electron-nucleus scat- 
tering has the form 



4^(1 + w e ^; 



1 + exp( 

where the constant c is such that the normalization / drn e (r) = A is fulfilled, 
and the parameters R e , a e , and the Fermi parameter W e can be found in 
Ref. [10]. For sufficiently heavy nuclei the radii are well described by the 
formula 



R e = (1.12A 1 / 3 - 0.86A- 1 / 3 ) fm, (2) 

while a e = 0.54 fm and W e = 0, which is the case of the heavy nuclei used at 
RHIC. From now on we set W e = 0, although the user may assign any value 
to this parameter (WF) in the input to GLISSANDO. 

In Glauber Monte-Carlo calculations one generates the positions of centers 
of nucleons in the nucleus of mass number A. A relevant point needs to be 
discussed here [11]. The nucleons are not point-like, hence they cause some 
washing out of the surface. Therefore the parameters R and a in the distribu- 
tion of centers of nucleons, n(r), must be such that when n(r) is folded with 
the charge distribution of the single nucleon, %(r), the distribution n e (r) is 
reproduced. We take the nucleon charge profile in the Gaussian form, 



where a = 0.79/y/S fm reproduces the nucleon rms charge radius of 0.79 fm. 
The folding gives 



n e (f) = J d 3 r' J d 3 r"n{f')n N {f")5 3 (f- f' - f"), (4) 
which for the distribution in the radial variable becomes 

n e (r) = 2ir J dz J r' 2 dr'n(r')nN(Vr' 2 + r 2 — 2rr'z). (5) 
Now, it turns out that the Woods-Saxon form of n e (r) is well reproduced with 
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a Woods-Saxon form of n(r), albeit with different parameters. A fairly good 
fit is achieved with the formula 

R = (1.12A 1/3 - 0.622A- 1/3 ) fm, (6) 
a = 0.46 fm, (distribution of centers). 

Note a shrinkage of the surface thickness parameter a and a slight increase 
of R, with an overall effect of a sharper distribution of the centers than the 
electron scattering profile. 

For the convenience of a user wishing to further study these issues, we include 
a Mathematica script prof ile_f olding.nb, where the folding of the Woods- 
Saxon profile with the nucleon profile is analyzed. This script is an additional 
tool and is not an integral part of GLISSANDO. 

There are two more effects that influence the nuclear density distribution in 
the Monte Carlo approach. The first one is the shift of the generated coordi- 
nates to the center-of-mass frame of the nucleus. This causes some shrinkage 
of the distribution. The other effect is the introduction of nucleon-nucleon ex- 
pulsion distance, which is a popular way to simulate the short-range repulsion 
in Glauber-like calculations. The centers of nucleons in each nucleus cannot 
be closer to each other than a certain expulsion distance d. The magnitude of 
d should be of the order of the hard-core repulsion range in the nuclear po- 
tential. Typically, values of a fraction of a fm are used. GLISSANDO offers two 
ways of generating locations of nucleons satisfying the expulsion constraint. In 
the first method, labeled "fix-last" , the positions are generated subsequently. 
When the currently generated nucleon is in the forbidden region, i. e. its cen- 
ter is closer than d to a center of any other previously generated nucleons, it 
is then generated anew. The procedure is repeated, until a "good" location 
is found. In the second method, labeled "return-to-beginning", when a nu- 
cleon is generated in the forbidden region the whole sequence of locations of 
the previously found nucleons is discarded and the procedure is started from 
the beginning, until a "good" sequence of nucleon locations is formed. The 
two methods yield virtually the same results for all investigated observables, 
and since the "fix-last" method is a few times faster (with d = 0.4 fm), it 
is preferred. Nevertheless, the user has a choice between the two generation 
methods. The value d = 0.4 fm is a typically range of the repulsive core in the 
NN interaction and is used as default in our input files. 

The repulsion implemented via the condition d > increases somewhat (at 
the level of 1%) the size R of the nucleus, which is noticeable in both above- 
described methods. This effect, and the previously mentioned shrinkage due to 
the shift to the center of mass, must be compensated by changing appropriately 
the parameters of the "bare" distribution from which the positions of centers of 
nucleons are generated. This is important, as otherwise a distribution different 
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from the desired one will be effectively obtained. We find that the following 
parameterizations do a fairly good job (here we give them for the "fixed-last" 
method only): 



R=(1.1UA 1/3 - 0.246A" 1/3 ) fm, 
a = 0.45 fm (d = fm — no expulsion) 



(7) 



R= (1.113A 1/3 - 0.277A" 1/3 ) 
a = 0.45fm (d = 0.4 fm), 



fm, 
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R= (1.103A 1/3 - 0.550A" 1/3 ) 



fm, 



(9) 



a = 0.455 fm (d = 0.8 fm). 

These formulas are implemented in GLISSANDO. In order to use them, one 
needs to set RWSA=-1 and RWSB=-1 in the input file. The parameterizations 
are provided only for the cases d = 0, 0.4, or 0.8 fm and for RET=0 (fixed-last 
method). Otherwise the user must explicitly give the values for the distribution 
parameters in the input file. 

We summarize: parameters flTHHj) produce approximately the Woods-Saxon 
distributions of centers of nucleons with parameters dSJ), which upon folding 
with the finite-size charge distribution of the nucleon reproduce to a very good 
approximation formula flTJ with parameters (j2J). 

The user can override the default parameters in the input file, as explained 
below, as well as examine the resulting distribution of centers of nucleons. 
The functional form of the distribution may be changed straightforwardly by 
modifying the formulas in functions . cpp. 

For the special case of the deuteron we use the Hulthen distribution [12] for 
the distance between the centers of the two nucleons, 



with a d = 0.228 fm" 1 and b d = 1.118 fm" 1 . 
2.2 Collision 

After generating the two nuclei, their centers-of- masses must be separated by 
the impact parameter b. Technically, the (x, y) coordinates of the nucleons in 



n d (r) 



1a d b d (a d + b d ) 
(a d - b d ) 2 



[exp(-2a d r) + exp(-2b d r) - 2 exp(-(a d + b d )r)] 



(10) 
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nucleus A are shifted by a fixed transverse vector such that their center-of- 
mass position is (bB/(A + B), 0), while the nucleons in nucleus B are shifted 
such that their center-of-mass is at (— bA/(A + B),0). In this frame we call 
the point (0, the geometric center of mass. Next, the wounded nucleons and 
the binary collisions are counted. 



2.2.1 Hard-sphere wounding 

In the hard-sphere wounding prescription, a nucleon from one nucleus is 
wounded if its center passes closer (in the transverse plane) to the center 
of any of the nucleons from the second nucleus than the wounding distance 



r 



V^tiA. ( n ) 



where a w is equal to the NN inelastic cross section at a given collision en- 
ergy. The standard implementation of the wounded nucleon model at the 
RHIC yj sjvjv = 200 GeV energy assumes that the inelastic cross-section in 
the nucleon-nucleon collision is [13] 

o~ w = 42 mb. (12) 

The cooresponding value for the LHC energies of ^/snn = 5500 GeV is a w ~ 
63 mb, while for the SPS energies of ^Jsnn = 19 GeV one uses a w ~ 32 mb. 

Similarly, a binary collision occurs if centers of two nucleons, the first one 
from the nucleus A and the second from B, pass closer to each other than the 
distance tq. This means that in this model binary collisions are counted as if 
occuring with the same inelastic cross-section as for wounded nucleons. 



2.2.2 Gaussian wounding 

Instead of the hard-sphere method, one may use a smooth function to deter- 
mine the probability distribution of wounding or a binary collision. An option 
in GLISSAND0 allows the user to use the Gaussian wounding function, 



p(r) = G exp (^-^Pj , (13) 

where r is the transverse distance between the centers of nucleons and the 
normalization is / °° 2irrp(r)dr = 7TTq = a w . The parameter G controls the 
central value of the profile. The default value is G = 0.92, which is taken from 
the experimental pp analysis of [14,15]. In order to switch on this feature the 
preprocessor option -D_gauss_=l must be used in the makefiles. 



8 



2.2.3 Displacement of sources 

Since the nucleons have a finite size, the location of the source may be some- 
what displaced relative to the center of the wounded nucleon or the mean of 
the centers of the two nucleons for the case of the binary collisions. GLISSANDO 
may include this effect by randomly displacing the x and y coordinates with 
a Gaussian distribution of width DWS for the wounded nucleons and DBIN for 
the binary collisions. Non-zero values of these parameters increase somewhat 
the size of the fireball, as well as reduce the eccentricity. This has also been 
noticed in later studies of Ref. [16]. Although the displacement effect is very 
physical, it is non-standard, hence is not included in the default parameters. 
If the user wishes to put it in, the typical values for the parameters DW and 
DBIN should be around 0.7 fm, a typical nucleon size. 



3 Models 

In this section we describe the models implemented in GLISSANDO. 

3. 1 The relative deposited strength 

The collisions between nucleons (wounding or binary collisions) result in de- 
position of a certain amount of the transverse energy (or entropy) at the 
location of the source in the transverse plane, which is then carried away by 
the produced particles. This quantity may be different for the wounded nucle- 
ons and the binary collisions, as well as can fluctuate from source to source. 
For the studies of fluctuations the absolute normalization of the strength of 
the sources is irrelevant, hence we use the measure proportional to the de- 
posited transverse energy, the relative deposited strength (RDS). For the 
wounded nucleon model RDS is just half of the number of wounded nucle- 
ons, which is a convention. For more involved models it is composed of the 
number of wounded nucleons and binary collisions, possibly with an overlayed 
statistical distribution of strength (see the following). 

3.2 Wounded nucleon model 

In the wounded-nucleon model the RDS of 1/2 is attributed to each point in 
the transverse plane at the center of a wounded nucleons. This complies to 
the original convention of Ref. [18], where the average multiplicity of produced 
particles is N pSuTt = N pp N w /2, with N pp denoting the average multiplicity in 
proton-proton collisions. 
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3.3 Binary collisions 



For binary collisions the RDS of 1 is attributed to each collision point, which is 
taken as the mean of the transverse coordinates of the two colliding nucleons. 



3.4 Mixed model 



A successful description of multiplicities at RHIC has been achieved with 
a mixed model [17], amending the wounded nucleon model [18] with some 
binary collisions [19,20]. In this wounded nucleon obtains the RDS 

of (1 — a)/2, and a binary collision the RDS of a. The total RDS averaged 
over events is then (1 — a)N w /2 + aNun- The fits to particle multiplicities of 
Ref. [20] give a = 0.145 for collisions at v /s^7 = 200 GeV, and a = 0.12 for 
^sj^ = 19.6 GeV. 

Thus the a (ALPHA) parameter in the code controls the proportion of wounded 
nucleons to binary collisions, with a = corresponding to the pure wounded- 
nucleon model, a = 1 to binary collisions only, and the intermediate values to 
the mixed model. 



3.5 The hot-spot model 

We also implement a model with hot spots in the spirit of Ref. [21], assuming 
that the cross section for a semi-hard binary collision producing a hot spot 
is small, o"bi n — 2 mb, which is a parameter of the model. Importantly, when 
such a rare collision occurs it produces on the average a large amount of 
transverse energy, with RDS equal to aa w /o- hin . Technically, the hot spots 
are implemented as follows: the binary collision is generated according to the 
criterion of Sect. 12.21 but is accepted with the probability o"bm/ow, i.e. the 
RDS of acr w /(Tbin is assigned to the hot spot at its position in the transverse 
plane. Otherwise the RDS is zero. All the wounded nucleons have the RDS 
assigned to (1 — a)/2. In this paper we label any model with a^ n < <r w the 
hot-spot model, although it should actually mean Obm <C cr w . The weight an 
an event is on average equal to (1 — a)N w /2 + aiVbm, the same as in the mixed 
model, but it can fluctuate considerably from event to event depending on 
how many hot spots are created. 
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Fig. 1. Histograms for the superposition distribution generated by GLISSAND0 from 
the Poisson distribution, k = 2, 20000 events. The figure is generated by running 
. /glissando_prof ile . exe input_prof ile_0 . dat and then using the wd branch in 
the density tree of the output file glissando . root. Horizontal axis: weight, vertical 
axis - probability distribution. Normalization arbitrary. Plot generated with ROOT. 



3. 6 Superposition 



Each source from the previously described models (wounded-nucleon, binary, 
mixed, or hot-spot) may deposit the transverse energy with a certain prob- 
ability distribution. To incorporate this effect, we superimpose a statistical 
distribution of weights over the distribution of sources. GLISSANDO has three 
options here, controlled by the parameter MODEL: - no superposition (all su- 
perimposed weights equal to one), 1 - superposition of the Poisson distribution, 
and 2 - superposition of the gamma distribution. 

The Poisson distribution generates the discrete weights according to the for- 
mula 



, K WK exp(-K) 12 3 

[wk)\ k k k 

This distribution has (w) = 1 and a(w) = 1/y/R,. 

The T distribution generates continuous weights according to the density 



w K 1 n K exp(— kw) r ^ . 

g{w,K) = — -f we 0,oo . (15) 

This distribution also gives (w) = 1 and a(w) = l/^fn. 
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The superposition distributions can be supplied independently for the wounded 
nucleons and binary collisions. The parameters k for both distributions are de- 
noted as Uw and Ubin, respectively. 

The preprocessor variable _weight_ controls the output of the superposition 
density to the output ROOT tree. The binary glissando_profile.exe writes 
out the superposition weights for the wounded nucleons. 

The fact that the chosen superposition distributions have (w) = 1 is at no loss 
of generality, since, as already stated, the individual RDS used for carrying 
the statistical averages can be normalized arbitrarily. The average RDS per 
event is still (1 — a) N w / 2 + aN hin and does not change due to the superposition 
procedure. 

For the general case we assign the RDS of (1 — a)w/2 for the wounded nu- 
cleons and au>o" w /o"bin for those binary collisions which are accepted with the 
probability cr hhl /a w . The variable w is generated with one of the distributions 
described above. The wounded-nucleon model corresponds to a = 0, the mixed 
model to a > and <7bi n = cr w , and the hot-spot model to a > and a^ n < <7 W . 



4 Fixed- and variable-axes quantities 

We use the standard convention for the axes of the reference frame: the z- 
axis is along the beam, the x-axis lies in the reaction plane, and the ?/-axis 
is perpendicular to the reaction plane. The azimuthal angle G [— n, ir] is 
measured relative to the y-axis, hence 

y = pcos<ft, x = psmcp, (16) 
where p = \Jx 2 + y 2 is the transverse radius. 

We refer to the analysis in the reference frame fixed by the reaction plane as 
fixed-axes (also called standard in the literature), and to the analysis where 
the particles in each event are translated to the center-of-mass frame and 
aligned with the major principal axis of the second harmonic moment of the 
fireball as variable-axes (also called participant). 

If the reaction plane is determined in each event (which of course can never 
be achieved precisely in the experiment [22,23]), one can choose the reference 
frame fixed by the reaction plane. The boost-invariant (in the central rapidity 
region) two-dimensional profile of the fireball density, f(p,<f)), is obtained by 
averaging the RDS distribution over the events belonging to a particular cen- 
trality class. Thus the normalization is / d<ppdpf(p, <fi) = average total RDS. 
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For instance, in the wounded nucleon model / d(ppdpf(p, 0) = N w /2. The sym- 
metry f(p,(f>) = f(p, 7r — 0), occurring also for unequal nuclei, excludes odd 
cosine and even sine components in the Fourier decomposition, hence in the 
general case of unequal nuclei we have 



f(p, 0) = f {p) + 2/ 2 (p) cos(20) + 2/ 4 (p) cos(40) + . . . 

+ 2 9l (p) sin(0) + 2g 3 (p) sin(30) + . . . (A^B). (17) 

For equal colliding nuclei the symmetry f(p, 0) = f(p, —(f)) eliminates the sine 
functions, 



f(p, 0) = f (p) + 2f 2 {p) cos(20) + 2/ 4 (p) cos(40) + . . . , (A — B) 



If the switch SHIFT is set to 0, the fixed-axes distance p = \fx 2 + y 2 is mea- 
sured from the center of the geometric overlap of the two nuclei, which is the 
traditional way. For equal nuclei it is just the mid-point between their centers 
of mass. If SHIFT=1, then p is computed in each event relative to the center 
of mass of the fireball, i.e., the center of mass of the RDS deposited in that 
event. The harmonic moments obtained from (I17lll8p . which we call fixed-axes, 
are also called "standard" in the literature. 

In order to have some convenient quantitative measures of the profiles of 
Eq. (ITS]) one introduces their radial moments 



J2npf l (p)p k dp 

€k > 1 J2npf (p)p*dp liyj 

The choice of the weighting power k is arbitrary, with the typical choice k = 2. 
Higher values of k would make the measure more sensitive to the outer region 
of the system, see Ref. [8] for a discussion of this point. We note that in the 
popular notation 



e s td = £2,2 = e. (20) 



Experimentally, the reaction plane cannot be determined accurately. Numer- 
ous methods have been developed to analyze the azimuthal asymmetry in 
heavy-ion collisions. The role of the purely statistical fluctuations caused by 
the finite number of particles has been described in [5-8,16,24-29]. As shown 
in these papers, the effects of the variable geometry in the initial stage of the 
collision lead to important effects, in particular to an increase of the shape 
eccentricity resulting in increased values of the elliptic flow coefficient t> 2 . 
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In the variable-axes calculations the coordinates of the elementary sources in 
the given event are first shifted to the center-of-mass of the fireball. One then 
computes the principal axes of the ellipse of inertia of the fireball generated in 
the event, which is twisted relative to the reaction plane coordinate system. 
The angle between the major half- axis of the ellipse and the y axis is given by 
the relation 



taW = 2 ( f - {x)i f v (21) 
var(y) — var(x) 

where var denotes the variance. The brackets denote the average over the 
sources in the given event. The angle 0* fluctuates from event to event. The 
superscript * indicates quantities averaged in such a way, that first in each 
event one computes from Eq. (l2Tj) the rotation angle 0*, then the rotation 
is performed to the current principal-axis system, and finally averaging over 
events is performed. The procedure results in the "variable-axes" density 



r(j>, 0) = n{ P ) + 2f*(p) cos(20 - 20*) + 2/ 4 *(p) cos(40 - 40*) + . . . 

+ 2gl(p) sin(0 - 0*) + 2g*(p) sin(30 - 30*) + . . . (A^B), (22) 

or for the case of equal nuclei 



f*{p, 0) = /o 0) + 2/ 2 * (p) cos(20 - 20*) + 2/ 4 *(p) cos(40 - 40*) + . . . 

(A = B). (23) 



We refer to the above density profiles as to the variable-axes profiles. 



In analogy to Eq. ffl9|) we introduce the variable-axes moments 



, j27T P r(p)p k dp 

k > 1 s2<K P rMp k d P ' [ } 

In the commonly-used notation for the variable-axes or participant deforma- 
tion parameter one has 



The profiles and moments for higher harmonics are suppressed, similarly to the 
fixed -cLXGS CcLSG. This is clear, as the higher harmonics are evaluated relative 
to the axes determined by maximizing the quadrupole moment. As a result, 
only a few moments are needed to effectively parameterize the profile. 
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The shift to the center-of-mass of the fireball results in a slightly more compact 
distribution /o(p) than fo(p) in the geometric center-of-mass frame, see Fig. [6] 
below. 

In Ref. [8] we show that many of the qualitative and quantitative features of 
the Fourier distributions displayed above result from purely statistical consid- 
erations. 

GLISSANDO generates and stores the fixed- and variable-axes moments, as well 
as the two-dimensional density profiles of the fireball and their harmonic de- 
composition. They can be later used in other studies involving the shape of 
the medium, such as jet quenching [8], setting the initial condition for cascade 
or hydrodynamic studies, etc. The included program interpolation shows a 
sample use of the profiles in a C++ code. 



5 Typical sequence of running 

5. 1 Installation 

It is necessary to have the ROOT package installed. After obtaining the GLISSANDO 
distribution the user should unpack 

tar -xzvf glissando_l_5.tar .gz 

and run the installation script 
sh install 

As a result, executable binaries are created. Five of them perform simulations: 
glissando . exe for the nucleus-nucleus (A+B) collisions, glissando_prot . exe 
for the proton-nucleus (p+A) collisions, glissando_deut . exe for the deuteron- 
nucleus (d+A) collisions, and also glissando_profile.exe, which in addi- 
tion generates the nuclear profile of nucleus A while running the (A+B) 
collisions. All these codes implement the hard-sphere wounding profile, see 
Sect. 12.2.11 The binary glissando_gauss.exe performs the (A+B) collisions 
with the Gaussian wounding profile, see Sect. 12.2.21 The wounding profile 
is controlled with the _gauss_ flag in the makefiles (-D_gauss_=0 for hard 
sphere, -D_gauss_=l for Gaussian), hence the user may easily implement 
his choice also in the codes for the (p+A) and (d+A) collisions. The code 
interpolation . exe provides an example of creating interpolating functions of 
the two-dimensional density profiles and their Fourier-decompositions, based 
on the stored histograms. The code retrieve.exe illustrates a use of the full 
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Table 1 

Typical output to the console from GLISSANDO. 

GLISSANDO ver. 1.6 ( \protect\vrule widthOpt\protect\href {http://arxiv.org.abs/0710.573lHl 
tested with ROOT ver. 5.08—5.18 

Start: Tue Feb 5 20:47:46 2008 



nucleus - nucleus collisions 



parameters reset in input file input.dat : 
EVENTS 100000 
ALPHA 
BMAX 20 

generates output file glissando .root 
random seed: 1202240866, number of events: 100000 
197+197, RA=6.43, aA=0.45, RB=6.43, aB=0.45, d=0.4fm 
wounded nucleon model : sig_w=42mb 
(binary collisions not counted) 
hard sphere NN collision profile 
window: b_min=0fm, b_max=20fm 
fix-last algorithm 

event: 100000 (100°/,) 

Quantities for the specified b, N_w, and RDS window: 

A+B cross section = 6436. 6mb, equiv. hard-sphere radius = 7.15687fm 

efficiency (accepted/all) = 51.2208% 

N_w = 103. 179+/-105. 16 

RDS = 51.5897+/-52.58 

eps_fix_2 (std.) = . 315127+/-0 . 316701 
eps_var_2 (part.) = . 475002+/-0 . 275256 
eps_fix_4 = 0.0872776+/-0. 281066 
eps_var_4 = . 240773+/-0 . 31374 

x_cm = -0.000206385+/-0. 482281, y_cm = . 00400375+/-0 . 732762 

Finish: Tue Feb 5 20:54:11 2008 
(Oh :6m :25s) 

event tree, optionally generated by GLISSANDO. Several useful ROOT scripts are 
described in Appendix O 

5.2 Running the simulations 

For simulations of the A+B collisions the running command has the syntax 
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Vglissando.exe [input_f ile] [output_f ile] 

and analogously for the p+A (glissando_prot . exe) and d+A (glissando_deut . exe) 
cases. When the input and output file-name arguments are absent, their de- 
fault values are 

input.dat - default input 
glissando . root - default output 

Typical input files are also provided with the distribution. The input param- 
eters and their defaults are described in Appendix [B] Thus we may simply 
type 

. /glissando . exe 

The output to the console shown in Table [T] contains basic information on 
the run. The subsequent lines give the info on the input: the version of the 
code, initial time, type of reaction, name of the input file used and the values 
of parameters reset from the default, the seed for the ROOT random-number 
generator, the number of events, the mass numbers of nuclei, the "bare" (see 
Sect. 12. ip Woods-Saxon parameters and the expulsion distance, the type and 
parameters of the model (wounded, binary, mixed, hot-spot), the window in 
the impact parameter (and the number of wounded nucleons and the value 
of total weight), the type of the algorithm for the short-range repulsion, the 
dispersion parameters for the location of sources, and the counter for events. 
The final output consists of the total nucleus-nucleus cross section in the given 
window, ctab, the equivalent hard-sphere radius defined as and the 

efficiency parameter, denoting the ratio of events where the nuclei collided to 
all the Monte-Carlo generated events. Next come averages of basic quantities 
with their standard deviations: the number of the wounded nucleons, binary 
collisions, hot spots, RDS, eccentricity e in the fixed- and variable-axes frame, 
€4 in the fixed- and variable-axes frame, the center-of-mass of the fireball. 
Finally, the execution time is printed. 

The user may first wish to verify the resulting nuclear density profile. For the 
case with no expulsion distance (d — fm) simply run 

. /glissando_prof ile . exe input_prof ile_0 . dat 

The file input_prof ile_0 . dat sets the bare (as explained in Sect. 12. If) Woods- 
Saxon parameters for the distribution of the nucleon centers to R = 6.44 fm, 
a = 0.45 fm (gold) and uses no short-range expulsion (d = fm). The code 
generates the ROOT output file glissando . root. The user should then enter 
the ROOT interpreter 

root 
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OLISSANDU ver. 1.5, 100000 events 
A=197 

Woods-Saxon fit to the radial distribution of centers of nucleons 



x10 3 




2 4 6 8 10 12 



OLISSANDU ver. 1.5, 100000 events 
A=197 

Woods-Saxon fit to the radial distribution of centers of nucleons 



x10 3 




2 4 6 8 10 12 



Fig. 2. The histogram of the radial density for the distribution of centers of nucleons 
in the gold nucleus for the case of no expulsion distance (top), and with d = 0.4 fm 
(bottom). Fixed-last algorithm, the bare parameters given in the text. 

and execute the script 

.x fitr.CC") 

The histogram of the nuclear density of nucleus A stored in glissando . root 
is fitted to the Woods-Saxon shape. The script generates the top panel of 
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GLISSANDOver. 1.5 

197+197, 1 events 
b=6.0fm 

mixed model with hotspots: a w =42.0 mb, a bj/ =1.0 mb, a=0.145 

overlayed Gamma distribution with parameters U =2.0 (wounded) and V . =2.0 (binary) 



fixed-axes density! | variable-axes density | 




Fig. 3. The spatial distribution of RDS generated in the single event in the hot-spot 
model. 

Fig. [2] and displays the parameters R and a, which are very close to the values 
of parameterization ([7]). 

For the case with expulsion distance of d = 0.4 fm the "bare" Woods-Saxon 
parameters are R = 6.43 fm, a = 0.45 fm. The running sequence is as follows: 

. /glissando_prof ile . exe input_prof ile_04 . dat 
root 

.x fitr.CO'") 

The resulting plot and the fit parameters are shown in the bottom panel of 

Fig. m 

In studies of reactions involving new nuclei the user may tune the values of R 
and a in such a way, that with a chosen value of d a proper density profile of 
nuclear centers is obtained. 

Now we are ready to run the simulations. First, let us look at a single event 
and make a snap-shot of the transverse-energy distribution in the x — y plane. 
We use here for illustration the hot-spot model and fix the seed for the Monte 
Carlo generator, such the result is repeatable. 

Vglissando.exe input_snap.dat 
root 

.x density. C("") 
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Table 2 

Centrality classes in the impact parameter 6, number of wounded nucleons N w , and 
RDS for the mixed model with a = 0.145, gold-gold collisions. 

b min: 0.0199137, b max: 19.3733 
N_w min: 2, N_w max: 391 
RDS min: 1, RDS max: 377.349 
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The produced plots are shown in Fig [3] (of course, for a single event the fixed- 
and variable-axes distributions are the same up to a rotation in the x — y 
plane) . 

Next comes a full-fledged minimum-bias simulation for the mixed model for 
Au+Au collisions with 500000 events: 

. /glissando . exe input _minbi as .dat 
root 

.x centrality. C("") 

The script centrality .C determines the centrality classes in the 6-parameter, 
in the number of wounded nucleons, N w , and RDS. Since the number of pro- 
duced particles is proportional to the RDS, the last case corresponds in essence 
to the determination of centrality via the multiplicity of produced particles. 
The output of centrality. C is generated to the console as well as to the file 
centrality.dat. The file contains the information on the centrality classes 
in steps of 5%, as shown in Tabled Graphical representation of the results is 
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GLISSANDO ver. 1.5 
197+197, 500000 events 
b=0.0- 22.0 fin 

displacement par.: DW=0.7fm, DBIN=0.7 fm 
mixed model: o u .= 42.0 nib, 0^=42.0 mb, a=0.145 




b[fm] 



GLISSANDO ver. 1.5 
197+197, 500000 events 
b=0.0-22.0fm 

displacement par.: DW=0. 7 fm, DBIi\=0. 7 fm 
mixed model: a u = 42.0 mb, <3 hji =42.0 mb, a=0.145 




GLISSANDO ver. 1.5 
197+197, 500000 events 
b=0.0- 22.0 fin 

displacement par.: DW=0. 7 fm, DBIN=0. 7 fm 
mixed model: a H .= 42.0 mb, 0^=42.0 mb, d=0.145 




RDS 



Fig. 4. Centrality as the function of the impact parameter b, the number of wounded 
nucleons N w , and the relative deposited strength (RDS). 
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shown in Fig. __ generated by the script. 

A major interest has been recently attracted by the fixed- and variable-axes 
eccentricity, as its event-by-event fluctuations are linked to the fluctuations of 
the elliptic flow coefficient t>2. The following scripts show the dependence of e, 
e*, and their scaled standard deviations Ae/e and Ae*/e* as functions of N w 
and b: 

root 

.x epsilon.CC") 
.x epsilon_b.CC'") 

The scripts ask for a rebinning parameter (a small natural number) which 
smooths out the plots. This is desired when the statistics is small. Originally, 
200 bins for the impact parameter are used. If the rebinning parameter is 
greater than 1, then this number of neighboring bins is grouped together. The 
result for the dependence on N w are shown in Fig. __ The horizontal solid line 
indicates the value \J^/^ — 1 ~ 0.52, which is the limit of Ae*/e* for most 
central events [8]. 

Similarly, the script dxdy . C shows the dependence of the event-by-event stan- 
dard deviation of the center-of-mass coordinates of the fireball on N w : 

root 

.x dxdy.CC") 

We may now select the centrality window of interest and perform a simulation 
for that centrality class. For instance, for c = — 5% in N w we run 

. /glissando . exe input_0_5.dat glissando_0_5 . root 

(note that for better book keeping we now use a name for the ROOT output 
file, which is stored for a potential later use). The choice of the centrality in 
N w needs a careful selection of the window in b. Of course, one may use a 
minimum-bias window such as b = — 25 fm, but this causes a waste of time, 
as many events are not recorded and the efficiency is low. For the present 
case, we have found from Table __ (or in centrality.dat) that the selected 
centrality window corresponds to N w = 327 — 394. Thus, we have set W0=328 
and Wl=1000 (of course, any number > 394 is good) in the file input_0_5 . dat. 
It turns out that in this case it is enough to use b = — 4.5 fm. In order to 
verify if this selection is proper, we may perform a test by executing 

root 

TBrowser a 
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GLISSANDO ver. 1.5 

197+197, 500000 events 
b=0.0 - 22.0 fm 

displacement par.: DW=0.7fm, DBIN=0.7fm 
mixed model: o H = 42.0 mb, <5 bUi =42.0 mb, a=0.145 



0.8 




GLISSANDO ver. 1.5 

197+197, 500000 events 
b=0.0 - 22.0 fm 

displacement par.: DW=0.7fm, DBIN=0.7 fm 
mixed model: o H = 42.0 mb, c biii =42.0 mb, a=0.145 




Fig. 5. Fixed- and variable-axes eccentricities and their scaled standard deviations 
plotted as functions of N w . The reb inning parameter was set to 1. 

Selecting glissando_0_5 . root and clicking on the leaf b generates a histogram 
of b. The supplied range in b should be sufficiently wide, such that the his- 
togram is not chopped off on the sides. 

Next, the scripts density . C and prof ile2 . C generate plots of the two-dimensional 
fixed- and variable-axes profiles in the x — y plane and their Fourier compo- 
nents. Note that large statistics is required for the higher Fourier profiles. 



23 



GLISSANDOver. 1.5 

197+197, 500000 events 
b=0.0 - 4.5 fm 

displacement par.: DW=0.7 fm, DB1N=0.7 fm 
N F =32S-1000 

mixed model: G w = 42.0 nib, G bjii =42.0 nib, a=0.145 




P[fm] 



GLISSANDOver. 1.5 

197+197, 500000 events 
b=0.0 - 4.5 fm 

displacement par.: DW=0.7fm, DBIN=0.7fm 
N K =328-1000 

mixed model: a t = 42.0 mb, a bjn =42.0 mb, a=0.145 




p[fm] 



Fig. 6. Fixed- and variable-axes Fourier profiles for c = — 5%. 

root 

. x density . C ( "glissando_0_5 . root " ) 
. x prof ile2 . C ( "glissando_0_5 . root " ) 

An example is shown in Fig. [BJ where we give the fixed- and variable-axes 
profiles for I = and 2 harmonics. 



24 



The two-dimensional fireball profiles and their Fourier components may be 
used in further analyses, for instance in the calculation of jet quenching, 
where we need the shape of the absorbing medium, in event-by-event hy- 
drodynamics, where the initial conditions undergo fluctuation, in studies of 
the elliptic flow, etc. For that purpose we also provide the template C++ code 
interpolation, exe which generates interpolated values of the functions from 
the stored two-dimensional or one-dimensional histograms. The typical run- 
ning is 

. /interpolation . exe glissando_0_5 . root 
Test of normalization 
Mean RDS= 311.818 

Integral over 2 Pi rho fO(rho) drho= 311.824 

Integral over (xyhist) f (rho_x,rho_y) drho_x drho_y= 311.869 

INTERPOLATION 

ver 1.0 

1 - One Dimensional Interpolation 

2 - Two Dimensional Interpolation 

Type 1 or 2 (any other - Exit) 
1 

ID Interpolation 

Give name of ID histogram 
cOrhp 

Give value of rho : 1.5 

Interpolated value at 1.5 equals 4.31756 

Once more?(y/n) 

n 

The available names of one-dimensional histograms are 

cOhp, cOrhp, c2hp, c2rhp, c4hp, c4rhp, c6hp, c6rhp, s3hp, s3rhp, 

where c indicates the cosine, s the sine moments (for unequal nuclei), the num- 
ber 0,2,... denotes the index of the Fourier harmonic, and r labels the variable- 
axes profiles. Similarly, the two-dimensional interpolation may be tested. The 
names of histograms for the fixed- and variable-axes two-dimensional densities 
are: 

xyhist, xyhistr. 

The code serves as a template example for interfacing the GLISSANDO results 
to other applications in C++. 
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GL1SSAND0 ver. 1.5 
197+197, 100000 events 
b=10.0- 18.0 fm 

displacement par.: DW=0.7 fm, DBlN=0.7fm 
mixed model: G H = 42.0 mb, a bji =42.0 mb, a=0.145 



| fixed-axes density~| | variable-axes density | 




GLISSANDO ver. 1.5 

197+197, 100000 events 
b=10.0-18.0fm 

displacement par.: D W=0. 7 fm, DBIN=0. 7 fm 
N„=3-ll 

mixed model: <J M = 42.0 mb, a biti =42.0 mb, a=0.145 



| fixed-axes density"! | variable-axes density | 




Fig. 7. Three-dimensional density profiles for centrality 80-95%. 

The analysis shown above may be repeated for other centralities, e.g. 

Vglissando.exe input_30_40.dat glissando_30_40 . root 

. /glissando . exe input_80_95.dat glissando_80_95 . root 

Figure [7] shows the two-dimensional profiles for c = 80 — 95% obtained with 
the command 

root 

.x density . C("glissando_80_95 .root") 
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We may carry out the calculation for other variants of the Glauber model, for 
instance for the hot-spot model with superimposed T distribution 

. /glissando . exe input_hotspot.dat hs.root 
root 

.x centrality.C ("hs.root") 
.x epsilon.CC "hs.root") 

The code can be run for different nuclei (for the case of light nuclei the pa- 
rameters of the nuclear distribution should be examined case by the 
global parameterizations do not work accurately in this case). 

. /glissando . exe input_S_Pb.dat 
root 

.x centrality.CC") 
.x epsilon.CC") 



In the case of different nuclei it makes sense to shift the fireball to its center 
of mass also in the fixed-axes case, which is made by setting SHIFT=1 in the 
input file. 

For the proton-nucleus collisions we execute 

. /glissando_prot . exe input_prot . dat 
root 

.x centrality.CC") 

and for the deuteron-nucleus collisions 

. /glissando_deut . exe input_deut . dat 
root 

.x centrality.CC") 

We can obtain predictions for the upcoming LHC experiment by increasing 
the a w cross section to 63 mb (which is the appropriate value for the en- 
ergy y/SNN = 5500 GeV [13]) and making an educated guess for the mixing 
parameter, a = 0.2: 

. /glissando . exe input_minbias_LHC . dat 
root 

.x centrality.CC") 
.x epsilon.CC") 
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Fig. 8. Scattered plot for the correlation of numbers of wounded nucleons and binary 
collisions for the minimum-bias gold- gold collisions at RHIC. 

6 Muliplicity fluctuations 

Histograms ntarg, nbinar, and nwei in the GLISSAND0 output file contain, 
respectively, the average number of wounded nucleons in nucleus B, average 
number of binary collisions, and average RDS, obtained for a fixed number 
of wounded nucleons in nucleus A. Similarly, ntarg2, nbinar2, and nwei2 
contain the scaled variances of these quantities. Such calculations are relevant 
for the multiplicity fluctuations at CERN SPS, where the number of wounded 
nucleons in nucleus A (projectile) is measured with the help of the VETO 
calorimeter. 



7 Other results 

Numerous other results may be obtained interactively by entering the ROOT 
interpreter and accessing the GLISSAND0 output file. As an example, we show 
in Fig. [8] the event-by-event correlation plot between the number of wounded 
nucleons and the number of binary collisions for the minimum-bias gold-gold 
collisions at RHIC (input file input_minbias.dat). The plot is generated by 
the following interactive sequence of instructions: 

root 

TBrowser a 

[select the output file from GLISSAND0] 
[select events tree] 
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[right-click on the events tree in the left window 

and select Start Viewer] 
[drop leaves nwAB and nbin on the x and y fields] 
[hit draw from the bottom menu of the Viewer] 



8 Using the full event tree 

Setting the input parameter FULL to 1 causes the generation of the full event 
tree, stored in the output GLISSANDO file in the full_event tree. A sample 
run can be made as follows: 

. /glissando . exe input_full.dat 

The saved information contains the x and y coordinates (X, Y) and the RDS 
(W) of the source, as well as the event number (KK). Technically, a structure is 
formed 

typedef struct {Float_t X,Y,W; Ulnt.t KK;} SOURCE; 
static SOURCE tSource; 

TTree *full_event; 

full_event = new TTree ( "full_event ", "full event"); 
full.event -> Branch("full_source" ,&tSource, "X: Y:W:KK") ; 

consecutively filled with information on each source, and written to the branch 
f ull_source. Other useful information, such as the impact parameter in the 
given event, is stored in the events tree. A template use of the full event tree 
is shown in the program retrieve.exe. 

. /retrieve . exe [glissando output file] 

The full event information may be processed "off-line" with other existing 
programs. The user should be warned that the FULL=1 option results in 
generating a lot of data, about 10MB/10000 events. Thus the option should 
only be used if the full results are to be used for analysis by another existing 
program. 



9 Summary 

We hope that with its flexibility and simplicity GLISSANDO will become a 
useful tool for the heavy-ion community. The fact that the code is publicly 
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available allows for check-ups, additions and improvements. We have provided 
examples of numerous applications: determination of the A+B cross-section 
and centrality classes, analysis of eccentricities of the fireball both in the fixed- 
and variable axes frames, study of event-by-event fluctuations or correlation 
of various quantities. 

The authors are grateful to Adam Trzupek, Janusz Krywult, and Grzegorz 
Stefanek for useful comments, to Constantin Loizides for helpful e-mail ex- 
changes concerning the fluctuations of the center-of-mass of the fireball, and 
to Adam Bzdak for a discussion of the Gaussian wounding profile. We also 
thank Vittorio Soma for carrying out tests on the MacBook computer. 
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A Content of the package 



file name 



description 



glissando.cpp 

functions. cpp 

functions, h 

interpolation. cpp 

retrieve, cpp 

install 

Makefile 

Make file, prot 

Makefile, deut 

Makefile .profile 

Makefile, inter 

Makefile. retr 

input.dat 

input_profile_0.dat 

input_profile_04.dat 

input_snap.dat 

input _minbias . dat 

input_0_5.dat 

input_30_40.dat 

input_80_95.dat 

input_S_Pb.dat 

input_prot.dat 

input_deut.dat 

input_hotspot.dat 

input_minbias_LHC.dat 

input_full.dat 

info.C 

fitr.C 

centrality.C 

epsilon.C 

epsilon_b.C 

dxdy.C 

density. C 

profile2.C 

label. C 

readme_run 

profile_folding.nb 



the GLISSANDO source file 

the function source file 

the functions header file 

template interpolation code source file 

template code for retrieving info from the full event tree 

the installation shell script 

makefile for glissando.exe, nucleus-nucleus collisions 
makefile for glissando_prot.exe, proton-nucleus collisions 
makefile for glissando_deut.exe, deuteron-nucleus collisions 
makefile for glissando_profile.exe, profile for nucleus A 
makefile for interpolation.exe 
makefile for retrieve . exe 
the generic input 

input for glissando_profile.exe, d = fm 

input for glissando_profile.exe, d = 0.4 fm 

input for generating a single event 

input for minimum-bias gold-gold collisions at RHIC 

input for centrality 0-5%, gold-gold at RHIC 

input for centrality 30-40%, gold-gold at RHIC 

input for centrality 80-95%, gold-gold at RHIC 

input for S-Pb collisions at SPS 

input for proton-nucleus collisions 

input for deuteron-nucleus collisions 

input for the hot-spot model 

input for LHC predictions for Pb-Pb collisions 

input generating the full event tree 

script giving information on the stored output file 

script for producing and fitting the nuclear profile 

script for centrality classes 

script for eccentricity vs. N w , etc. 

script for eccentricity vs. b, etc. 

script for center-of-mass coordinates vs. N w 

script for two-dimensional profiles 

script for Fourier profiles 

script generating the label used by other scripts 
short instructions for typical running 
Mathematica file for folding the profile 
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B Description of input and output 



The basic model parameters, collected in Table IB. II can be supplied in the 
input file. The sign # at the beginning of the line comments out the line and 
then the default value of the parameter is used. See the sample file input . dat. 

The variables stored in the output GLISSANDO file are explained in Tables IB~2l 
and El 



Table B.l 

Parameters of the input file. 



name 


default 


description 


ISEED 





seed for the random number generator, if a random seed is generated 


EVENTS 


50000 


number of generated events 


NBIN 


40 


number of bins for histograms in p, x, or y 


FBIN 


72 


number of bins for histograms in the azimuthal angle 


NUMA 


197 


mass number of nucleus A 


NUMB 


197 


mass number of nucleus B 


RWSA 


6.43 


Woods-Saxon radius for the distribution of centers, nucleus A [fm] 






gold with the fix-last method 


AWSA 


0.45 


Woods-Saxon width, nucleus A [fm] 


RWSB 


6.43 


Woods-Saxon radius for the distribution of centers, nucleus B [fm] 


AWSB 


0.45 


Woods-Saxon width, nucleus B [fm] 


WFA 





the w parameter of the Fermi distribution, nucleus A 


WFB 





the w parameter of the Fermi distribution, nucleus B 


CD 


0.4 


closest allowed distance between centers of nucleons [fm] 


SNN 


42. 


AW "wounding" cross section [mb] 


SBIN 


42. 


AW binary cross section [mb] 


ALPHA 





- wounded, 1 - binary, 0.145 - PHOBOS for RHIC@200 GeV 


MODEL 





- constant superimposed weight=l, 1 - Poisson, 2 - Gamma 


Uw 


2 


Poisson or Gamma parameter for wounded 


T T1 " 

Ubm 


2 


Poisson or Gamma parameter tor binary 


DW 


0.7 


dispersion of the location of the source for wounded nucleons [fm] 


DBIN 


0.7 


dispersion of the location of the source for binary collisions [fm] 


WMIN 


2 


minimum number of wounded nucleons to record event 


WO 


2 


minimum allowed number of wounded nucleons 


Wl 


1000 


maximum allowed number of wounded nucleons 


RDSO 





minimum allowed RDS 


RDS1 


100000 


maximum allowed RDS 


GA 


0.92 


central value of the Gaussian wounding profile 


SHIFT 





1 - shift the coordinates of the fireball to the cm. in the fixed-axes case 






- do not shift 


RET 





- fix-last algorithm 






1 - return-to-beginning algorithm for nuclear density 


FULL 





1 - provide the full information on events (coordinates and RDS of sources) 






- do not 


DOBIN 





1 - compute the binary collisions also for the case ALPHA=0 






- do not compute the binary collisions for the case ALPHA=0 


BMIN 


0. 


minimum impact parameter [fm] 


BMAX 


25. 


maximum impact parameter [fm] 


BTOT 




range parameter for histograms [fm] 
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Table B.2 

Histograms stored in the output ROOT file. (.) denotes the mean and var the variance 
of the specified quantity. 



xyhist 
xyhistr 

cOhist 

cOrhist 

cOhp 

c2hp 

c4hp 

c6hp 

cOrhp 

c2rhp 

c4rhp 

c6rhp 

s3hp 

s3rhp 

nx 

nx2 

ny 

ny2 

neps 

neps2 

nepsp 

nepsp2 

nuni 

nepsb 

neps2b 

nepspb 

nepsp2b 

nunib 

ntarg 

ntarg2 

nbinar 

nbinar2 

nwei 

nwei2 



fixed-axes density in the x — y variables 
variable-axes density in the x — y variables 

fixed-axes density in the p — 4> variables (not normalized) 
variable-axes density in the p — 4> variables (not normalized) 
fo 

/2 - fixed-axes cosine harmonics 

h 

h 

fS 
/! 
ft 
fS 

93 
9*3 



variable-axes cosine harmonics 



fixed-axes sine harmonic 
variable-axes sine harmonic 
(x) [fm] vs. N w 
var(x) [fm] 2 vs. N w 
(y) [fm] vs. N w 
var(y) [fm] 2 vs. N w 
(e) vs. N w 
var(e)/(e) 2 vs. N w 
(e*) vs. N w 
var(e*)/(e*) 2 vs. N w 
event multiplicity vs. N w 
(e) vs. b 

var(e)/(e) vs. b 
(e*) vs. b 

var(e*)/(e*) 2 vs. b 
event multiplicity vs. b 
(N*) vs. 

varW)/(JV|) vs. N£ 
(N hiQ ) vs. 

va,i(N hin )/(N Un ) vs. N£ 
(RDS) vs. N£ 
va,v(RDS)/(RDS) vs. 



C Description of the ROOT scripts 

Each script is called with the name of the ROOT file generated by GLISSANDO, 
root 

.x <script>.C(" [{\tt GLISSANDO} file]")} 

The empty string in the argument of <script>.C("") is equivalent to the 
default <script> . C ( "glissando . root " ) . 

info.CC [output ROOT file]") 
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Table B.3 

Trees and their contents stored in the output ROOT file. 



TTree param all parameters of the calculation 



TTree density (generated only by glissando_profile.exe) 

r radius of the nucleon in nucleus A [fm] 

wd weight generated by the superposition distribution 



TTree phys 

sitot the nucleus-nucleus cross section [mb] 

eps_fixed event-by-event average e 

eps_variable event-by-event average e* 

sigma_eps_fixed event-by-event standard deviation of e 
sigma_eps_variable event-by-event standard deviation of e* 

TTree events 

nwA number of wounded nucleons in A 

nwB number of wounded nucleons in B 

nwAB total number of wounded nucleons 

nbin number of binary collisions 

npa RDS 

b impact parameter 

es fixed-axes epsilon (standard), < r 2 cos(20) > 

ep variable-axes epsilon (participant), < r 2 cos(2(0 — 0*)) > 

es3 < r 2 cos (30) > 

es4 < r 2 cos (40) > 

ep4 < r 2 cos(4(0 - 0*)) > 

es6 < r 2 cos (60) > 

ep6 < r 2 cos(6(0 - 0*)) > 

es3s < r 2 sin(30) > 

phir the rotation angle 0* 

xx x cm. coordinate [fm] 

yy y cm. coordinate [fm] 



produces info on the stored ROOT file. 
fitr.CC" [output ROOT file]") 

produces the density profile of centers of nucleons and its Woods-Saxon fit. 
The script works on ROOT files generated with glissando_prof ile . exe, which 
contain the treel tree an the leaf r. 

central ity.CC [output ROOT file]") 

produces centrality windows for the stored ROOT file. This makes sense for 
minimum-bias calculations. 

epsilon. C(" [output ROOT file]") 
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produces plots of the mean and scaled standard deviation of the deforma- 
tion parameters e and e* as functions of N w . The script asks for a rebinning 
parameter (a natural number). 

epsilon_b.C(" [output ROOT file]") 

same as epsilon.C but with N w replaced by the impact parameter b. 
dxdy . C (" [output ROOT file]") 

produces the event-by-event standard deviation of the coordinates of the 
center-of-mass of the fireball as a function of N w . 

density. C(" [output ROOT file]") 

produces plots of the two-dimensional fixed-axes and variable-axes spatial dis- 
tribution of RDS. 

prof ile2.C(" [output ROOT file]") 

produces plots of the fixed-axes and variable axes Fourier profiles of spatial 
distribution of RDS. 
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